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We investigate the role of nonperturbative quark-gluon dynamics in soft high energy 
processes. In order to reproduce differential and total cross sections for elastic proton- 
proton and proton-antiproton-scattering at high energy and small momentum transfer 
it turns out that we need two scales, the gluonic correlation length and a confinement 
scale. We find a small gluonic correlation length, a ~ 0.2 fm, in accordance with recent 
lattice QCD results. 
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Introduction. Soft hadronic collisions at high energy are described successfully using Regge 
and Pomeron phenomenology Jl], 0. Yet an outstanding problem is the more basic under- 
standing of pomeron exchange in terms of the underlying quark-gluon-dynamics. 

At large center of mass energies a/s but small momentum transfers \t\ ~ (1 — 2) GeV 2 , one 
expects that non-perturbative features of QCD, such as gluonic vacuum properties, enter as 
important ingredients in hadron-hadron scattering amplitudes. In this context several recent 
papers || HI pj have emphasized specifically the role of gluonic field strength correlations in 
hadronic scattering processes at high energy. 

The aim of the present paper is to follow up on this theme and gain further insight into 
the basic scales which govern soft high energy collisions. In particular, the question arises 
whether the small gluonic correlation length recently found in lattice QCD calculations || 
is compatible with the scales deduced from the t-dependence of differential cross sections 
for high energy proton-proton and proton-antiproton collisions. Our starting point is the 
assumption that the scattering of composite hadrons can be reduced to the scattering of 
their quark constituents; we develop a model of quark-quark scattering in terms of non- 
perturbative gluon dynamics. Our basic picture is that soft interactions between quarks 
at high energy are described by the propagation of projectile quarks in a nonperturbative 
gluonic environment produced by the target quarks. The quark sources themselves are tied 
together by confining forces inside the hadron projectile and target. 

Propagation of High Energy Quarks in a Gluonic Background. We start from the 
QCD Lagrangian 

C(x) = -hrG^(x)G^(x) + $(x) [i^D^x) - m ] ${x) (1) 

with the usual field strength tensor G^ = d fl G u — d u G^ + ig[G^, G u ] and the covariant 
derivative D^(x) = d^ + igG^x). Here G^(x) = G®(x)T a are the gluon fields contracted with 
the SU(3) co \ or generators T a = X a /2. The quark field i[)(x) includes three color components, 
and mo is the current quark mass. Consider the vacuum-vacuum transition amplitude 

Z = JVGV^V^ exp ^ J d*x C(x)} , (2) 

where the integration is over all independent quark and gluon fields. Our aim is to reduce 
the complexity of the gluon field configurations appearing in eq. (2) to the point where one 
can see the explicit appearance of gluonic correlation functions. 

Let f[G, tp] be a functional of the gluon and quark fields. We define an average over the 
gluon fields (G-average) as 

(f[G,j>]) G := J VG exp J d 4 x [-^trG^(x)G^(a;)]}/[G,^. (3) 

For a pure Yang-Mills theory this is simply the expectation value of f[G}. With the help of 
eq. (3) the amplitude (2) can be rewritten: 
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Z = J Dip Dip exp |i J d 4 xip(x) [ry^c^ — m ] ip(x)j ■ 
■(V exp{-^|d 4 x Ja M (^)^(x)}) G , 

(4) 

where j£(x) = ip(x)j^T a ijj(x) is the quark color current. Here V denotes Dyson pathordering0. 
In the next step we use the cumulant expansion to reexpress the G-average in (4) as 
follows: 



G 



'V exp|-z 5 j<&cf a Gl 

00 (—i a ) m r r 
ex P{ E —zr~ / d ^ ■ ■ ■ / d^m^(n) • • .Jt(%) < G£(*i) • • • G£(s m ) >J. 

m=l m - J J 

(5) 



The index c now includes Dyson-ordered cumulants. For a given function /, these are 
constructed out of the G-averaged moments by collecting the connected parts: 



< f(x) > c = < f(x) > G , 
< f{x 1 )f(x 2 ) > c = < f(x 1 )f{x 2 ) > G - < f(xi) > G < f{x 2 ) > G ,etc. 

(6) 

Eq. (5) is exact for commuting fields and correct for noncommuting gluon fields if one 
includes moments up to second order. For higher orders eq. (5) is still a good approximation 
if the length scale of the underlying correlations is small M compared to the size of the 
interacting composite particles. 

We abbreviate the quark-gluon interaction part on the right hand side of the cumulant 
expansion (5) by — ig J d^cj^(x) < G*(i) >, with the effective gluonic background field 



[G;(x)) =: <G*(*)> G + 
1i (m + l)\ 



E T^f^j / dV " 7 d W£(*0 • • -CW < Gl{x) . . . G£ (x m ) > c 

(7) 



The corresponding effective Lagrangian 

£ eff (x) = ^(ar) [i-fdpix) - m ] ^(z) - <?j£(a;) < G°(x) >, (8) 

describes the propagation of a quark in the presence of a gluonic background created by 
other quark sources. 
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In the following we restrict ourselves to lowest non- vanishing order in the expansion (7), 
which gives for the interaction part of eq. (§): 

4ff (*) = / ( G » Gb M) c ■ ( 9 ) 

Note that this approximate expression still incorporates important features of non-perturbative 
gluon-dynamics. Its effects are twofold: first, it dresses single quarks, shifting the current 
quark mass m to (average) constituent masses m. At high energies this effect is negligeable 
because all masses are small compared to the center of mass energy. Secondly, it generates 
interactions between the color currents of projectile and target quarks. 

Quark-Quark-Scattering. In the next step we construct the quark-quark scattering am- 
plitude in a high energy approximation. Using standard LSZ reduction we write down the 
corresponding 4-point Green function and arrive at the scattering matrix: 

(PzPa\ S - 1 |pip 2 ) = i (2vr) 4 5 4 (p 4 + p 3 - p 2 -pi)M qq 

= jd% /dV. exp {i („ .*. + „. P3 )} etoxtow & - m )(i ft. - ro )*(*„ 

(10) 

The full two particle wave function ^{x 3: x±;pi,p2) is defined as: 

»(*,,*;*.,») = /dV/A, exp { - i (, 2 . P2 + x 1 .p 1 )}(0|^(x 4 )«,3)«x 2 )«x 1 )|0> 

(i ^ 2 +m)M(p 2 )(« ^ +m)u(pi). (11) 

The free quark wave functions are <pp(x) = u(p) exp {— ip ■ x}, with Dirac spinors u(p). 
At high energy, we are able to factorize ^(x 3 , x 4 ; pi, p 2 ) into a product of single particle wave 
functions, <p pi (x 3 ) and ^(x^), each moving in an external field created by the other particle, 
provided that we assume the color matrix commutator [T a , to vanish. We have convinced 
ourselves that corrections due to the non-commuting part amount to at most a few percent 
for the leading interaction term. With the approximation (^|) this external potential can be 
constructed explicitly. 

To arrive at the factorization we introduce lightcone momentum variables p± = p° ± p 3 . 
We neglect backscattering and make use of the fact that at large momenta p the only non- 
vanishing current matrix elements are < p± \j± \p± > while < p T \j± \p T > and the transverse 
components vanish. The momenta of the scattered quarks differ from p + or p_ only by a 
small transverse amount |p_i_| <s: \p\. Couplings to sea quarks are suppressed in this limit. 
We can now solve the Dirac equation for the scattered quarks in the high energy approxima- 
tion of ref. Q and arrive at the following eikonal expressions for the quark wave functions 
(j) pi (x 3 ) and P2 (x 4 ): 

<f) pi {x) = exp|i^da/ + V_(x / + ,aj±,a:_)|0j x (a;) =: <$>-{x)<f) pi {x), (12) 
(j) P2 (x) = exp|z£^dx'„V + (x + ,a ;i .,a; / _)|0° 2 (x)=:$ + (x)</)° 2 (a;). (13) 
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Here the components of x are the lightcone variables x+, x±, X-. The quark interacts with the 
effective gluonic background field along its lightlike path x + (or and we have introduced 
the operators and $ + which distort the quark wave functions along their trajectories. 
The gluonic background is described by the potentials V±(x): 



V-(x) = -*^T a /dV^ a (y)7-^(y)(G^(*)Gt(y)) c , (14) 



2 

V + (x) = -ijT a j d%^ 1 (y) 1+ T b ^ Pl (y)(G a + (x)G b _(y)) c . (15) 

Each of these potentials incorporates the full wave function of the second quark. Therefore 
the coupled equations flT2|) and (^) have to be solved self-consistently. 
Introducing cm. coordinates l/2(x3+£4) and relative coordinates z = £3 — £4 in eq. QTO) one 
obtains the scattering amplitude M. qq . At high energy we approximate j9~ (7+cL + 7_<9 + ). 
We can then perform the integration along the relative coordinates z + and z_, and we are 
left with an integral over the transverse coordinates b = z±. The resulting quark-quark 
scattering amplitude is 

M gq(<l 2 ) = ^u(p 3 )j + u(p 1 )u{p i )j^u(p 2 ) J d 2 6 exp{-iq-b} 
tr [$_(oo, 6/2, 0) - 1] tr [$+(0, -6/2, 00) - 1] , 



(16) 



where a color trace has been taken. 



Gluonic Correlation Function. The potentials V±(x) involve the gluonic correlation 
function < G±{x)G T (y) > c . With the help of the nonabelian Stokes theorem [[K] 



we can 



transform the integrations over the paths along the gluon fields into surface integrals over 
the field strengths G a : 



S exp|-^T a T fe J da +± J da' 1 (Gt ± (x)G b +± (y))^ ^°\x) 



(17) 



Here S refers to surface ordering which is the remnant of the former pathordering, da ±± 
are the surface elements. To identify our correlator with the gauge invariant correlation 
function < g 2 G'] lu (x)(j)(x, y; xa)G b (y) > c we have to work in the Fock-Schwinger-gauge, 
where the Schwinger string <ft(x, y; xq) = 1. At this point we loose manifest local SU(3)-color 
gauge invariance and keep only global SU(3)-color symmetries. Note that there is also a 
dependence on the reference point xq. This dependence can be shifted systematically into 
the higher order cumulants |§. 

The correlation function \ G a ^ v {x)G b pa {y)\ is our main input. Using Lorentz and translational 
invariance we can write it as ITTII 



G%(x)G b pa (y)) r = (G« u (x)G b ^(x)) r [g, p g ua - g^g vp ) D f ^-Ji ) (i 8 ) 



a 
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where D(\x — y\/a) is a function of the distance of the correlated fields. The most general 
expression includes also a derivative term which fulfills the abelian Bianchi identities, whereas 
the nonabelian character is described by the nonderivative part ||. In lattice calculations 
the nonabelian term (|18|) is found to dominate 0. For the nonlocality D(\x — y\)/a) we use 
a parametrization given by Dosch et al. |L3| 



„f\x-y\\ f d 4 k f i, , A I 27tt 4 k 2 \ 

D K*) = iw ^ h k ■ (' - »)} (-— ( g- 9^/64)' ) ' (19) 

normalized to D(0) = 1. This introduces the gluonic correlation length a which enters in 
the calculation of cross sections. The diagonal elements (g 2 G a ^ v {x)G^ u {x) ) are related to 



the gluon condensate | 12| . The nondiagonal elements vanish when taking the color trace. 



Note that in leading order the potentials V±, eqn. 15), are proportional to the gluon 
condensate times the non- locality D which scales like a 4 . 

Differential Cross Sections. Given the quark-quark scattering matrix A4 qq , we still need 
to know the transverse distributions of projectile and target quarks. We describe those by a 
properly normalized form factor F(q 2 ). We introduce the eikonal phase, or profile function 

x{b2) = I exp {iq ■ b} F V)AW<? 2 ) (20) 

and construct the scattering matrix for the composite hadrons: 

M(q 2 ) = 2 J d 2 b exp {-iq ■ b} [ exp {ix(b 2 )} - l] . (21) 
The differential cross section is 

and the total cross section satisfies the optical theorem 

a tot = lmM(t = 0). (23) 

For the hadronic formfactor we choose a monopole form F(t) = A 2 / (A 2 — t) with A = v^6/^ C) 
where r c plays the role of a confinement radius. Taking a dipole or a gaussian formfactor 
instead would not change the value of r c very much but would result in a slightly less optimal 
fit to the cross section at — t > 2 GeV 2 . 

In table 1 we list our best-fit parameters, which reproduce total and differential cross sections 
for pp and pp scattering. We find good agreement with the ISR pp data MM as well as with 



the pp data at ^fs = 546GeV Hl5 



Summary and discussion 

We have presented a description of high energy elastic hadron scattering using a model 
characterized by two basic scales: the gluonic correlation length and a confinement radius. 
The gluon correlation length determines the interaction of constituent quarks in the projectile 
with those in the target. The strength of this interaction is governed by the gluon condensate 
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for which we find values remarkably close to results obtained in the QCD sum rule approach 
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The gluon correlation length a turns out to be primarily responsible for the dip in the 
differential cross section. Its value a ~ 0.2 fm is small in comparison with typical hadron 
sizes and agrees with recent SU(3) lattice results on the gluon field strength correlator |J. 
It is also compatible with the small size of the scalar glueball wave function extracted from 
lattice calculations [HJ. We note that the dip in da/dt shows up in the amplitude already 
before eikonalization. 

A confinement scale r c ~ 0.5 — 0.6 fm completes the input to obtain a good overall fit to 
the measured differential cross sections for both pp and pp. Both length scales, a and r c , 
together determine the slope of da/dt at small \t\. This slope is frequently interpreted in 
terms of effective mean square radii of the interacting hadrons: 

b = (J^ ln ~^J = \( <r2 > P r °j ectile + < r<1 > targot) • (24) 
Phenomenological analysis |17| emphasizes the apparent equality of such hadronic radii with 



the empirical electromagnetic radii. The present study points to a different interpretation. 
Dosch et al. [Tj| , in a similar investigation, give larger values for both the correlation length 
a and the confinement radius r c . For the latter they take the value of the proton charge 
radius, r c = 0.86 fm, and find a = 0.35 fm. We can reproduce their results by calculating the 
slope parameter b using our scattering amplitude, but restricted to the Born approximation, 
i.e. the lowest order contribution only. We point out that the full calculation, with the basic 
quark-quark interaction iterated to all orders, induces changes of roughly 30% in comparison 
with the Born amplitude which result in the smaller values of a and r c given in table 1. 





figure la 


figure lb 


figure 2 




pp 


pp 


pp 




52.8GeV 


62.5GeV 


546GeV 


a (fm) 


0.196 


0.199 


0.220 


r c (fm) 


0.54 


0.54 


0.53 


< (a/n)G 2 > (GeV 4 ) 


0.014 


0.014 


0.008 



Table 1: Gluonic correlation length a, confinement scale r c and gluon condensate < (a/7r)G 2 > 
used to reproduce the data for pp scattering at yfs = 52.8 GeV and 62.5 GeV and pp scattering 
at y^s = 546 GeV as shown in figs. 1,2. (The "standard" value of the gluon condensate is 
< (a/vr)G 2 >= 0.012 ± 0.003 GeV 4 ||/.J 
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Figure 1. pp-scattering at \fs = 52.8 GeV (a) and 62.5 GeV (b). The lines present our 
optimal fits with the parameters given in table 1. Fig. la shows in addition a fit with 
a = 0.3 fm.The data are taken from 
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Figure 2. pp-scattering at = 546 GeV. The line presents our optimal fit with the param- 
eters given in table 1. The data are taken from [15 ). 
(Appended as uu-encoded files) 
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Figure 1. pp-scattering at yfs = 52.8 GeV (a) and 62.5 GeV (b). The lines present our 
optimal fits with the parameters given in table 1. Fig. la shows in addition a ht with 
a = 0.3 fm.The dat taken from [14]. 
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Figure 2. pp-scattering at y^i = 546 GeV. The line presents our optimal fit with the para- 
meters given in table 1. The dat taken from [15]. 



